Stability of Active Transmission Lines 
with Arbitrary Imperfections 

By H. E. ROWE 

(Manuscript received August 23, 1963) 

Two sufficient conditions for the stability of one-dimensional active 
transmission lines with arbitrary imperfections (i.e., discrete or continuous 
reflections) are derived. The first stability condition guarantees stability 
for any arbitrary distribution of reflection. The second stability condition is 
restricted to a special case of interest that includes discrete reflectors with 
nominally equal magnitude and spacing; the stability condition for this re- 
stricted class is greatly improved over the general stability condition de- 
scribed above. 

These results, aside from their own interest, provide rigorous justifica- 
tion for previous calculations for the gain statistics of such a device with 
random discrete reflectors. 1 They may also be used to find an upper bound 
on the probability of instability of such a device with random reflectors. 

Certain types of optical maser amplifiers and traveling-wave tubes pro- 
vide examples of practical devices with distributed gain to which these re- 
sults, or similar ones, might be applied. 

I. INTRODUCTION 

The preceding paper 1 has considered the theory of active transmis- 
sion lines with discrete imperfections. First, lines with equally-spaced 
identical reflectors were studied; in particular, gain-frequency curves 
were determined as functions of the various parameters, and the sta- 
bility of the device was studied under these special conditions. It was 
pointed out that the mathematical expression for gain would yield a 
perfectly definite result for any values of the parameters, but that this 
mathematical result would have physical significance only if the device 
is stable, i.e., does not oscillate. 

Next, the case of random imperfections was studied. 1 Here the statis- 
tics of the transmission were determined in terms of the statistics of the 
discrete reflectors, which were assumed to have random position and 
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magnitude. Again, these results have physical significance only if the de- 
vice is stable (or if the probability of instability is negligible) . However, 
in the random case no precise information about stability was given; 
the computed statistics of the transmission were felt to be valid if the 
rms magnitude of the discrete reflectors was sufficiently small, but only 
intuitive feelings of what was "small enough" were available. 

In the present paper we derive a sufficient condition for stability of an 
active transmission line with arbitrary reflectors; we further show (by 
one example) that this sufficient condition cannot be greatly improved 
(if at all) in the general case. This result gives useful information re- 
garding the range of validity of the calculations of the preceding paper 1 
for the transmission statistics of active transmission lines with random 
reflectors. This general bound on stability may be improved if additional 
information is known about the distribution of reflectors; one such case 

of interest is treated. 

The mathematical model chosen for this problem is discussed in detail 
in Ref. 1. A line with N discrete reflectors is shown in Fig. 1 (which is 
identical to Fig. 1 of Ref. 1). The wave traveling to the right at distance 
z is denoted by Wo(z), the wave traveling to the left by Wi(z) ; W (U-+ ) 
and Wi(L k +) are the right- and left-traveling waves just to the right 
of the fcth reflector, as indicated in this figure, while W (L k -) and 
Wi(Li : -) are the right- and left-traveling waves just to the left of the 
feth reflector. 

In the absence of reflections the forward and backward waves vary as 



Wo(z) « e 

Wi(z) « e 



— forward wave 

— backward wave 
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Fig. 1 — Line with N discrete reflectors. 
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where 



T = -a + jp, a > 0. 



(2) 



The line has gain, so that a > 0. From (12) of Ref. 1, the wave matrix 
for the cascade connection of the fcth line section of length h and the 
kth reflector is 



X k = 
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f" e +n " 
_+jc k e~ 


+rifc-| 




Vl - c k 2 
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= X,- 


pr (L*+)~| 

_TF x (L,+ )_ 



cJ ^ 1, 



(3) 



(4) 



where | c* 1 is the magnitude of the reflection coefficient for the fcth re- 
flector. The over-all transmission matrix for the entire line of Fig. 1, 
denoted by X, is given by the matrix product of (13) of Ref. 1 : 



i = n^, 



'w (oY 

Tfi(0) 



= X- 



'W (L N +Y 
Wi(Lg+) 



(5) 



(6) 



For convenience, denote the elements of the over-all transmission 
matrix X as in (14 ) of Ref. 1 . 



X = 



.l'21 X 2 2 



(7) 



X is given by (3) and (5). Assume the device is operated as an amplifier 
with matched input and output; setting Wi{L N -\-) = 0, the complex 
transmission gain G T is given by 



Gt = 



W»(L N +) 
PTo(O) 






(8) 



Now .Tn is a function of T and of all of the k's and c*'s. We may concep- 
tually investigate stability in the following way. Imagine that c* is re- 
placed by eCk throughout this analysis; e is a variable parameter that 
scales the magnitudes of all of the coupling coefficients. Let e be increased 
from 0, and for eacli value of e examine :v n [which in (8) is the reciprocal 
of the transmission gain, and so may be regarded as the transmission 
loss] as a function of frequency co (or of the phase constant /3, which is 
assumed proportional to frequency, since the line is distortionless) 1 over 
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the entire range — °° < w < +«. We determine in this way the mini- 
mum value of | x n | for each value of e. As e increases, this minimum value 
of | x n | will eventually just drop to zero, for a critical value of e which 
we denote by e c . Thus, as e -» e c the gain | G T | -» °° for a particular 
value of co, and the device oscillates. e c is the dividing line between sta- 
bility and instability; if e c > 1, the original device, with the parameters 
Ck and h , is stable. 

Such calculations have actually been carried out in Ref. 1 for devices 
with identical, equally-spaced reflectors. In this case the gain G T is a 
periodic function of frequency co, so that only a finite portion of the fre- 
quency axis (i.e., one period) must be investigated. In general, however, 
G r is not periodic; since we cannot investigate numerically the entire 
co-axis, it is not obvious how to investigate stability for the general case. 

In the remainder of this paper we determine a sufficient condition 
that guarantees the stability of a general active line with arbitrary dis- 
crete imperfections. In particular, consider such a device, illustrated in 
Fig. 1, characterized by (3), (5), and (6), with arbitrary a, c* , and l k . 
We show below that any such device satisfying the condition 

£ tanlT 1 \a\ < 2 sinir' — — (9) 

i^i V2 

must be stable. Many practical devices will have large gain, and hence 
must have small reflections. In such cases e~ aLli « 1 and | c t - 1 « 1 ; under 
these conditions a slightly poorer stability condition derived from (9) 
is useful. 



n r P - aL N 

Y,\ci\ < tanh 2 sinir 1 —— 

fci L V2 



do) 



In the high-gain case the right-hand side of (10) may be made simpler 
still by further degrading this stability condition. We may show, for 
example, that 

tanh [2 sinh- 1 e -^\ ^ 0.932 V2e~ aLN , 8.686«L, V ^ 10 db. (11) 
V2 J 

Thus a slightly poorer version of (10) is 

£ I a I ^ 0.932 -s/2 e~ aL \ 8.686 aL N ^ 10 db. (12) 

The stability condition of (12) is valid when the one-way gain of the 
active medium exceeds 10 db. As the lower bound on the one-way gain 
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of the active medium increases beyond 10 db, the numerical factor 0.932 
on the right-hand side of ( 12) increases, approaching 1 as the lower bound 
on the gain approaches infinity. This is readily seen from (10); as 
aL N — > oo t e~ aLN — > 0, so that the sinh -1 and tanh functions in (10) may 
be approximately replaced by their arguments for sufficiently large aL N . 
However, direct calculation with (10) is straightforward; the result of 
(12) (or similar equations) is intended principally to illustrate the gen- 
eral behavior. 

Thus (9) or the successively poorer versions of (10) and (12) guaran- 
tee that the device will be stable, even for the worst possible choice of 
the c k and l k . Equations (9), (10), and (12) are each sufficient, but 
not necessary, conditions for stability. These results are derived in 
Sections II, III, and IV. In addition, a better bound is obtained for a 
special case in which the reflection coefficient is distributed more or less 
uniformly with distance z along the active line, in a certain sense to be 
described more precisely in Section V below; these results include many 
cases of interest. Finally, some numerical examples illustrating the use 
of these two different types of bounds are given in Section VI. 

II. DIFFERENTIAL EQUATIONS EQUIVALENT TO MATRIX RELATIONS 

Consider the following differential equations: 

W '(z) = -rir ( 2 ) + jr(z)W 1 (z), 

WY(*) = -jr(z)W (z) + TW.iz). 

These relations have the form of the coupled line equations with a gen- 
eral continuous coupling coefficient. In the present case, W (z) and TT 7 i(z) 
are the right- and left-directed traveling-wave complex amplitudes, and 
r(z) is the continuous reflection that couples the two waves to each other. 
Equation 13 is readily obtained as a limiting form of the matrix relations 
of (3), (5), and (6) by assuming very small, closely spaced discrete re- 
flectors whose magnitude varies slowly with distance. Thus in the matrix 
relations of Section I above set 

h = Az. (14J 

Assume that c k varies slowly with k. Then we set 

c k = r(kAz)-Az, (15) 

where r(z) is a continuous function. We now let Az —* so that the num- 
ber of discrete reflectors — * <*> ; during this process the continuous func- 
tion r(z) is fixed and the c k determined by (15), so that the magnitudes 
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of the individual reflectors — > as Az — » 0, Then the matrix relations of 
(3), (5), and (6) will yield the continuous differential equations of (13). 
The analysis is straightforward and quite similar to that of Ref. 2 for a 
similar problem, and so will not be given here. The above discussion of 
(13) as an appropriate limiting continuous form of the matrix relations 
of Section I is given only to provide some physical motivation for con- 
sidering (13), and plays no part in the mathematical analysis to follow. 
The case of isolated, discrete reflectors, characterized by (3), (5), and 
(6), may conversely be regarded as a special case of continuous reflection 
in (13), in which the continuous reflection r(z) becomes a sum of suit- 
able 5-functions, one located at each discrete reflector. Thus we show 
that if r( z) in (13) is given by 



r( z ) = Xtanh l d-8(z - Li) 



(16) 



where in Fig. 1 L, is the total distance from the input of the line to the 
ith reflector, then the solutions to (13) at the output of the line, i.e., 
W (L N +) and Wi(L N +), are given in terms of the input conditions 
IFo(O) and Wi(0) by (3), (5), and (6). 

Consider the typical fcth section of line, of length h , followed by the 
fcth discrete reflector, as illustrated in Fig. 1. In the line section between 
the (fc - l)th and the fcth reflectors r(z) = 0, from (16). Therefore in 
this region the solution to (13) has the form of (1); the forward and 
backward waves are uncoupled, and have the same propagation constant. 
We may thus write the solution between the (k — 1 )th and kth reflectors 
in the matrix form 



[ 



jf (l*-i+) 

TFi(L fc _H-)J 



r e+ n k 











'W (L k -Y 
WiiLt-) 



(17) 



where W(L k — ) indicates a wave amplitude evaluated just to the left of 
the fcth reflector, W(L k +) just to the right. 

We next evaluate the transmission matrix for the fcth reflector, i.e., 
the fcth S-fimction of (16). This calculation may be performed by setting 



L k < z < L k + A. 



tanh 1 Ck 
r(z) = " _ ^ _ ' 

0, otherwise. 



We then determine the matrix T(A), 
~Wo(L k + A) 



_W 1 (L k + A)_ 



= T(A) 



W (L k ) 
DFi(I*)J 



(18) 



(19) 
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Then as A -> 0, r(z) -* tanh ! c ir 8(z - L k ), and lim T(A) = T(0) 

A-0 

yields a matrix relating the wave amplitudes W and W y on the two sides 
of the fcth S-function of r(z) [see (16)]. This analysis is again similar in 
motivation, although different in detail, to that of Ref. 2 for a similar 
problem. Since r(z) in ( 18) is constant throughout the region of interest, 
( 13) becomes a linear differential equation with constant coefficients, and 
is readily solved by the usual techniques. The solution for general A may 
be written in matrix form, yielding 2'(A) of (19), as follows: 



T(A) = 
1 


f-K-e r ^~ + K + e~ rAV ~ 


e v 

rr TAV — 

K+e v 


-e~ r * 


K + - K. 


- K-e 



-rA,/- 



K± = -j 



, 1 ± V 



tanh -1 c k 



K+K- = 1 



TA 

tanh -1 



c k 



1 



J 



TA 



K+ - K- 2 V 

v~ = 



V 1 + (t)' 



Taking the limit as A -^ 0, (20) -(23) yield 



w (L k +y 

Wt(L k +) 



= r(0) 



Wi(I*-) 



where 



T(0) = lim T(A) = l 

Vl - c? _-j Ck 



JCh 

1 



Inverting (24), 



where, from (25) 



w (L k -) 

W x (L k -) 



= T~ l (0) 



T~\0) = 



VT^ 



w (L k +y 

W x {L k +)_ 

r i -jcki 

c * 2 _+jc k 1 



(20) 

(21) 

(22) 

(23) 

(24) 
(25) 
(26) 
(27) 
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From (17), (26), and (27) we now have 



'Wo(L k -i+y 



- x k - 



w<>(L h +y 



(28) 



where X k is as given in (3). Equation (28) is identical to (4). Finally, 
the solution to (13), with r{z) given by (16), is given by (3), (5), 

and (6). 

The equivalence of (13) and (16) with (3), (5) and (6) is useful be- 
cause the original matrix problem may thus be regarded as a special 
case of a pair of differential equations. Stability appears to be more 
readily studied for the more general continuous case described by the 
differential equations; these results may then be applied to the special 
discrete case of interest here. 

III. SOLUTION BY SUCCESSIVE APPROXIMATIONS (PICARD's METHOD) 

We summarize the solution of (13) by successive approximation, fol- 
lowing the same general approach as in Ref . 3 for a similar problem. First, 
it is convenient to make the following transformations: 

Wi(z) = e +rz -C n (z). 
Substituting (29) into (13), we have 

G' (z) = jr(z) e^'GM 
(&'(«) = -jr{z) e- 2r 'G (z). 

Assume that the device is operated as an amplifier with matched input 
and output. It proves convenient in the following analysis to take the 
input at the right-hand end of the amplifier, i.e., at z = L N , where L N 
is the total length, and the output at the left-hand end, i.e., z = 0; this 
is just opposite to the choice made in Ref. 1 and in Section I above 
[particularly in (8)]. The useful output is then the left-directed traveling 
wave at z = 0, i.e., Wi(0), corresponding to an input taken to be the 
left-directed traveling wave at z = L N , Wi(L N ). Since the device is 
matched at both ends, W (0) = 0; W (L N ) ^ 0, since this quantity 
corresponds to the reflected wave at the input end (i.e., at z = L s ) 
of the amplifier. 

Now assume for convenience a unit-amplitude output wave: 

Wi(0) = 1. (31) 
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As noted above, since the output is matched, 

JFo(O) = 0. (32) 

We seek Wi(Lk), the input corresponding to the output of (31); since 
unit output has been assumed in (31), the complex transmission gain 
G, will be 

Gt = ur it \ ' (33) 

where Wi(L N ) is the solution to (13) subject to the initial conditions of 
(31) and (32). 

The transmission gain is readily stated in terms of the solutions to 
(30), which were obtained from (13) via the transformation of (29). 
Thus, consider (30) subject to the initial conditions 

C'o(0) = 0, 

(34) 
Gi(0) = 1, 

obtained from (31) and (32) via (29). The complex transmission gain 
Gy of the amplifier is then given by 

Gr =e- lL ».-±— (35) 

(j i{L/ N ) 

where (n{L N ) is the solution to (30) subject to the initial conditions of 
(34). 

We now seek the solution to (30), with the initial conditions of (34), 
via Picard's method of successive approximations. 4,5 Assume the 
(n — l)th approximation to the solution is available; let us denote this 
approximation by 6'o (n _i)(2) and G X{n _i ) (z). Then the (n — l)th 
approximation is substituted into the right-hand side of (30) and the 
right-hand side integrated to yield the nth approximation. 

Go(»)(z) =j I r(s) 6 +2r 'Gi( n _i)(s) ds. 

(36) 
G Hll) (z) = 1 -j / r(«)c- fr 'G (,-„(«)d«. 

We take the initial (0th) approximation as simply the initial conditions 

of (34): 

G om (z) = 0, 

(37) 

^'i(0)(z) = 1. 
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Writing 

GoooCs) — G (n-i)(z) = 0000(2), 
G Un) {z) - Gi(„-n(z) = 0100(2), 

we have 



(38) 



(39) 



Gooo (2) = Jldomiz), 
k=i 

G Un) (z) = 1 + J2guk)(z). 

k -1 

From (36) and (38), the 0's of (39) are given as follows: 

floooU) - 3 f r(s) •**W»(f) ds, n ^ 1. (40) 

■'o 

giooiz) = -j I* r(s) e - 2r Vo (n -i)(s) ds, n * 1. (41) 

Jo 

gom(z) = 0, giM = 1. (42) 

From (40) -(42) 

0000(2) = 0, n even. 

0100(2) =0, n odd. 

Thus only odd terms appear in the top summation of (39), and only 
even terms appear in the bottom summation of (39). 

We next obtain bounds on the magnitudes of the terms in the series 
of (39), thus showing that these series converge as n — > » for all finite 
z, so that the solutions to (30) subject to the initial conditions of (34) 
are 

-r. 

G (z) = Z 0000(2), 
Gi(z) = £ 0100(2), 

n=0 

with 0000(2) and 0100(2) as given by (40)-(42). The analysis is sug- 
gested by that of Ref. 3. We show that: 



(43) 



= 0, n even. 

00(h)(2) I 

< L J ° ' 

n\ 



< V± J_ , n odd. 
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[j[ | ,(,) I d,]' 



n even. 



9un)(z) 



< e 



n\ 



= 0, n odd. 

where from (2) 

T = -a + j(l, a = -Re Y > 0. 
Suppose that (40) is true for some even n. Then from (40) 

,« \ f \r(s) Ids'] 
| go <n+1) (z) | £ / | r(t) | e*« e+ 2at Ul =L 

Jn M.I 



303 



(46) 



(47) 



<fl 



■sr[( ,rw| 4{/. ,|r(i)| *] (48) 

[[\r(s)\ d s] n+1 



(n+l)! 
in agreement with (45). Substituting this result into (41), 

r \f\r(s)\ds] ' 

Jo (n + 1)! 

+2« r ,= 2 p .< -|„+l 



tfl 



_ (49) 



■I / |r(s) Us 



= e 



+2az 



j[ I r(«) | ds] 



rc + 2)! 



in agreement with (46). Noting (42) and (43), the results of (45) and 
(46) hold for all n by induction. 

The bounds of (45) and (46) guarantee the convergence of the series 
solutions of (44) under quite general conditions. It is readily seen that 



| Go(z) | ^ sinh / | r(s) | ds 



■]■ 



I Gi(z) | ^ e +iaz cosh ["jT | r(s) | <fel . 



(50) 
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The series solutions of (44) converge for all finite z, so long as the con- 
tinuous reflection coefficient is absolutely integrable, 

f \r(s) \ds < oo. (51) 

Jo 

In particular, note that r{z) may contain 5-functions, as in (16),, so that 
the above bounds may be applied directly to the discrete case of Sec- 
tion I. 

The solutions to (30) given by (44) and (40) -(43) thus converge for 
all finite z in the case of interest. However these formal mathematical 
solutions have physical significance only when the device to which they 
apply is stable, i.e., does not oscillate. In the following section we use the 
bounds of (45) and (46) to obtain a sufficient condition guaranteeing 
stability in the general case. 

IV. BOUNDS ON STABILITY — GENERAL CASE 

Consider a general amplifier described by (13) or equivalently by 
(30). Assume the total length is given by L N . We may investigate 
stability as indicated following (8). Replace the continuous reflection 
coefficient r(z) by e-r(z), where e is a numerical parameter. Let e be 
increased from 0, and for each value of e determine the maximum value 
of the transmission gain | G r | as a function of frequency o>. From (35) 
the maximum value of \G T \ corresponds to the minimum value of 
| Gi(Ltf) |. As e approaches a critical value, denoted above by e c , 
| Gr |max -» °° and | G^L*) | mi „ -> 0; if e r > 1 the original device is 
stable. 

From (40)-(44), 

00 

Gi(L N ) = 1 + E gun)(L N ). (52) 

n=2 
n oven 

Noting that r(z) has been temporarily replaced by e-r(z), for sufficiently 
small e a lower bound on the magnitude of Gi{L N ) is given by 

\Qi{L K )\ ^ 1 ~ E |0ioo(WI • (53) 



-2 

n even 



Both sides of (52) and (53) are functions of frequency w, through their 
dependence on the propagation constant /3. Using the result of (46) in 
(53), 



Gx{Lti) | fc 1 - E ' 

n=2 
7i even 



r L N "1" 

2aL N UP 'J (54) 
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Since the expression on the right-hand side of (54) is independent of 
the propagation constant /3 and hence of the frequency w, this expression 
is also a lower hound on | Gi(L s ) |min , the minimum value of | Gi(L N ) \ 
as a function of w. 



Gi{Ls) | min £ 1 - S e 



2ah N 



\ e-r(s) ds] 

.J* =L (55) 



As e increases from 0, the lower bound on | G\{L N ) |min given by (55) 
steadily decreases, and for some particular value of e ^ e c approaches 0. 
Therefore if 



E 



fl'N ~|n 

I | € r(s) |rfsj 



< e 



-2a L N 



(56) 



stability is guaranteed. If (56) is satisfied for e = 1, then stability is 
guaranteed for the original amplifier, with reflection coefficient r(z). 

Consequently, a sufficient stability condition for an active transmis- 
sion line with a general continuous reflection coefficient r(z), described 
by either (13) or (30), assuming the device to be matched at both ends, 
is given by 



E .2 



'..v 



r(s) \ds 



1 



n 



< e 



2a L N 



(57! 



This may be written 

r f LN 1 

cosh / J r(s) I ds — 1 < e 



-2a L N 



or further 



sinh 2 



/ I r(s) I ds 

In 



2 



< 1 .•*- 



Finally, taking the square root of both sides of (59) we obtain 



sinh 



I r(s) I ds 



< 






(58) 



(59) 



(60) 
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or eqiiivalently 



L 



\r{s) \ds < 2 sinh -1 —7=- (61) 



as sufficient conditions for stability for a general active transmission line 
with an arbitrary continuous reflection coefficient r(z). 

We may now apply the result of (61) to the discrete case of Section I 
above by making use of the results of Section II. As noted in Section II, 
if the continuous coupling coefficient r(z) is a series of 5-f unctions of 
the form given in (16), then the solution to (13) is identical to that for 
the discrete case, given in (3), (5), and (6). Since the stability condition 
of (61) holds true in general, it may be applied to the discrete case by 
substituting (16) into (61), yielding 

2 tanbT 1 | a | < 2 sinh -1 ^^ . (62) 

Equation (62) is a sufficient condition for stability for a general active 
transmission line with arbitrary discrete reflectors, having reflection co- 
efficients Ci located at arbitrary positions along the line. Equation (62) 
is the result stated in Section I as (9). This inequality is a sufficient con- 
dition for stability; if the inequality is satisfied, the device must be 
stable. This condition is not necessary for stability; many devices that 
violate (62) or (9) are stable. 

The weaker bounds of (10) and (12) are readily obtained from the 
basic result of (62) or (9) by straightforward use of inequalities. From 
(62) or (9) we must have 

tanh -1 \d\ < 2 sinh -1 ^^ i - 1, 2, • • • N. (63) 

Since the function y = tanh -1 x is concave upward for x > 0, 

tanh -\ T < tanh X _*m <x <x m <l. (64) 

Xm 

Therefore, from (63), 

. e~ aL ff 
2 sinh —7= 

tanh" 1 I Ci I < f —ZTTn 'I « I • (65) 



re~ aL N 
2 sinh- 1 —^ 
V2j 



tanh 2 sinh -1 
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Therefore if the relation 

AL T e~ aL N~\ 

J2\ci\ ^ tanh 2 sinh" 1 — y^ (66) 

is satisfied, then the condition of (62) must also be satisfied, so that (66) 
is a slightly poorer sufficient condition for stability; this result was given 
in (10). Finally, since the function y = tanh [2 sinh -1 x] is concave 
downward for x > 0, 



i , _ tanh [2 sinh x n 



tanh [2 sinh x] ^ - --x, g x £ x m . (67) 

X m 

As a particular instance let us choose x m = (l/\/20) — 0.2236; then 
(67) becomes 

tanh [2 sinh -1 x] ^ 1.863 x, ^ x ^ —7= - 0.2236. (68) 

By using (68) to decrease the right-hand side of (66), we obtain the 
slightly poorer sufficient condition for stability 






(69) 
- 0.932 -s/2 e~ aLN , 20 logio e"' N ^ 10 db 



given in (12). 



V. BOI'NDS ON STABILITY — SPECIAL CASE, INCLUDING REFLECTORS OF 
NOMINALLY EQUAL MAGNITUDE AND SPACING 

The bounds on stability derived in Section IV in the general case 
guarantee stability for the worst possible arrangement of reflectors. 
Thus in many cases the sum of the magnitudes of the reflectors may far 
exceed the bound given by (9), (10), 01 (12) without causing instability. 

These general bounds guarantee stability even if we have no informa- 
tion whatever about the distribution of reflectors. If we do have such 
additional information, it should be possible to make use of it to find im- 
proved bounds. As a trivial example, in the treatment of equally spaced, 
identical reflectors in the previous paper 1 exact stability conditions were 
obtained; we will see in Section VI that for this case the sum of the 
magnitudes of the reflectors at the boundary of instability may far ex- 
ceed that given by (9), (10), or (12). 

In the present section we consider a somewhat restricted special case 
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in which the reflection coefficient is almost uniformly distributed in a 
certain sense. We assume that 

R-(z -/) £ f \r{s) \ds ^ R-(z + g), 

J ° (70) 

R > 0, f^O, g ^ 0, 

where R, /, and g are constants. Equation (70) states that the indefinite 
integral of the absolute magnitude of the reflection coefficient is con- 
strained to lie between two straight lines of the same slope R, separated 
by the horizontal distance h given by 

h=f + g, h^0. (71) 

It turns out that the final bounds of this section are better the smaller 
the separation h. This is to be expected, since the smaller the separation 
of the two straight lines given by the right- and left-hand sides of (70), 
the more constrained is the reflection coefficient r{z). 

The presence of sufficient length of perfect (i.e., reflectionless) active 
line at either end will needlessly increase /and hence h in (70) and (71), 
and hence needlessly degrade the final stability condition given below. 
Such a length of perfect line cannot affect the stability, but merely alters 
the gain of the device (assuming it is stable). Therefore for purposes of 
the present stability analysis sufficient lengths of perfect active line 
should be removed from each end so that h is minimized, and hence the 
best possible bound is obtained. Removal of any additional lengths of 
perfect active line from either end will do neither good nor harm to the 
final stability condition. 

A few examples serve to illustrate the general nature of the restriction 
of (70). First suppose that r(z) is equal to a (positive) constant, 

r(» = r . (72) 

Then (70) is true with 

R = r 

/ = 0, g = (73) 

h a / + g - 0. 

The separation h [of (71)] between the straight lines of the two sides of 
the inequality of (70) is zero in this case. Equations (13) or (30) are 
readily solved exactly for the reflection coefficient of (72) by slight modi- 
fication of the results of (18)-(23), in particular by first replacing 
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tank -1 Ck — * r A and subsequently replacing any remaining A's by 
A — > L, where L is the total length, in these equations. From this exact 
solution precise stability conditions may be obtained for the case of 
constant (continuous) reflection coefficient; we expect the bounds of the 
present section to agree with this exact result when we set f = g = 0. 

Similarly, the parameters of (73) apply to the bounds of (70) when 
the (continuous) reflection coefficient is a square wave of constant 
absolute value r , with arbitrary transitions between the +r and the 
— ro sections. 

The above two examples utilize a continuous reflection coefficient. 
However, our particular present interest lies in some of the discrete cases 
of the preceding paper. First, consider the case of identical, equally- 
spaced reflectors of Section II, Rcf. 1 ; the relations of (70) are illustrated 
for this case in Fig. 2. A less-restricted case is provided by the case of 
reflectors of identical magnitude but random spacing, where the fluctua- 
tion in spacing is very small compared to the average spacing, treated in 
Section III of Ref. 1. The relations of (70) for this case are shown in 
Fig. 3 ; the randomness in spacing has resulted in a slightly wider separa- 
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Fig. 2 — Identical, equally spaced reflectors. 
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Fig. 3 — Identical, randomly spaced reflectors. 

fcion than in Fig. 2 between the dashed lines that enclose the staircase 
curve of 



f\r(s) 



ds. 



Since in this case the magnitudes of the reflectors are strictly constant, 
the "risers" of the staircase have the same size, while the "treads" vary 
in length. It is clear that if the magnitudes as well as the spacings of the 
reflectors vary slightly, both the "risers" and the "treads" of the stair- 
case will vary slightly, but otherwise the behavior will be much the same 
as in Fig. 3, so that the restriction of (70) may be satisfied with small 
separation between the straight-line bounds. 

While the discrete cases of the preceding paragraph, which have re- 
flectors of nominally equal magnitude and spacing, are of principal in- 
terest here and supply the motivation for the analysis of the present 
section, discrete reflectors having quite different distributions from the 
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above may also fall within the restriction of (70) with small separation 
of the bounding lines; one such case is illustrated in Fig. 4. (Note that 
reflectors of both signs are indicated in the lower drawing of this figure, 
by S-f unctions with both positive and negative magnitudes.) 

The above cases, which satisfy the restriction of (70), may be re- 
garded as having the absolute magnitude of the reflection coefficient 
more or less constant in a certain sense, in that 



f \r(s)\ 



ds 



is approximately proportional to z [see (70)]. Thus we seek bounds on 
stability in the case of ( 70) that are similar to those obtained for constant 
reflection coefficient [see (72)]. 

We again use the solution by successive approximation given in Sec- 
tion III above. The discussion of (29)-(43) remains appropriate for our 




r(z) 



t(z) = sum OF 6- FUNCTIONS 



Fig. 4 — More general case satisfying the restrictions of (70). 
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present purposes. However, greatly improved hounds over those obtained 
in (44)-(51) may be obtained because of the additional restriction of 
(70) imposed in the present section; in contrast, the bounds of (44)-(47) 
of Section III hold true in general, and specifically when the restriction 
of (70) is not satisfied. 

Consider the series solutions of (44). From (42) 



g m (z) = 1, 0o(o>(z) = 0. 
Note also (43). We show that: 

loooGOl < R2 Q a + h ) 2e2< " 

[ *a+»)Ms-)'i 

n even, n ^ 2. 



(74) 



(n/2)-l 



gi(nM I = 0, 



n odd. 



(75) 



I 0oc> («) I = 0, 
|*o>(f)| <R Q^ + h ) 



fe#+M*l 



n even. 



(n-l)/2 



B 1 ^ 



(^ 



n odd. 



In (75) and (76), R and h are the parameters of (70) and (71). 
First, from (40), (42) or (74), and (47), 

| 0b(i)(f) | ^ / o " <T 2< " | r(s) | rfs = JT' e*" d jf' | r(«) | d*J 

j" I r(0 | A + 2a J' e- 2a> |"|' | r(0 | <ftl ds, 



= e 



(76) 



(77) 



where we have made use of integration by parts. Using (70) in (77), 
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gom(z) I ^ e 2a '-R(z + g) + 2aR f e~ 2as ■ (s + g) ds 
- e~ 2ai -R{z + (,) 



[ 2a 



ze- 2c " + g(l-e 



■■»] 



(78) 



= i (1 " '"*"> + R <> < Ya + R(J " + g) > 

where in the final step we have used the fact that/ ^ 0. Finally, substitut- 
ing the definition of h from (71) into (78), 



0o(n (z) | </2^ + *). 



(79) 



Equation (79) agrees with (76) for n = 1. 
Next, from (41), (47), and (79), 

I <7i(2) (z) I < R (~ + hj f e +2ae | r(s) | ds 



1 

2a 



= R [ ^ + h ) I e +2as d 



r(t) 



I*] 



*(£ + *]' 



(80) 



•' \r(t)\dt 

Using (70), (80) becomes 

\9un(*)l <R 2 Q^ + hje 2a! -(z + g) 

- R 2 (± + hj 2a f\ 2as (s - /) ds 

= R*(± t + h)e i "-(z + g) (81) 
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Finally from (71), (81) becomes 

\9wW\ <#(& + *)' *~> (82) 

which agrees with (75) for n = 2. 

We now establish the bounds of (75) and (76) by induction. Suppose 
that (75) is true for some even n ^ 2. Then from (40) and (47), 



gou>+i)(z) | < 



■(£+0 



(n/2)+l 



R n hf- 



(I- 



(83) 



II 



where 



'-£[• + 6- O'l^lX 1 * 01 *]- (84) 

Integrating (84) by parts, 

-6-)n-6-)'r 

•N |r(0|<ft|<b- 

Using (70) and (71), we have from (85) 

l(n/2)-l 



(85) 



/^ [^ + (|- iW R(z + g] 



6-0f [• + G- 1 )*J" r '" <- - y »* 



-«/.' («-/) <*[«+(!- 1)»] 



(n/2)-l 
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(80) 



= Rft 2 



+ 



fl 



#*s 



£ - 1 ft 



B 



(1) 



GO 



[■ + G-o»r+H' + G-o*n- 



where the last step follows from the preceding one because n ^ 2 [from 
(75)],/ ^ [from (70)], and h ^ [from (71)]. Using the inequality 



x k + «B* l < [x + (e//c)] fc , z ^ and e > 0, (87) 



(86) yields 
/ 



© L ^ ; J (i) 



[z + w/i] n/2 . 



(88) 



Substituting (88) into (83), 



/,'■ [± + h) (* + nA)l 



n/2 



<7o<» + i>(«) | < 22 (^ + A] ! = ^ -^- J- . (89) 



(t> 



Recalling that // is some even integer ^ 2 in (89), (89 ) agrees with (76). 
Next, from (41) and (47), using the result of (89) 



<7i(»+-2)(z) I < 



R " + ii + >) 



(n/2)+l 



J 



n 
2f 



(90) 
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where 



(91) 



(92) 



(Is 



J = f e +2a '(s + nh) n,2 d \J | r(t) \ dt\. 

Integrating (91) by parts, 

J = e^iz + nhy' 2 ^ \r(t)\dtj 

- 2a f e +2aB (s + nh) nl * [£ |K0 I dt\ ds 

-\l e +2o *(s + nh) {nm - 1 [jf I r(t) | <ft] ds. 

Using (70) and (71), we have from (92) 
J ^ e 2a \z + nh) n,2 R(z + g) 

- R2a f e +2a8 (s + nh)" 12 (s -f) ds 
Jo 

\2> Jo 

= e laz (z + nh) n,2 R(z + </) 

- R f (s-f)d [e +2aB (s + nfc) n/2 ] 
•'o 

= e 2w (2 + nh) nl2 R{z + g) - R(z - f)e 2a! (z + nh) nl2 

- Rf(nh)" ri +r[' e +2a °(s + nh)" 12 ds (93) 

Jo 

= Rhe 2a: (z + nh) n/2 - Rf(nh)" 12 

+ * f(, + nfc) B/s d( e to ) 

Za Jo 
= Rhe- a \z + nh)"' 2 - fl/(«A)"" + £ e 2a '(z + n/i) n/2 

_ * (nh )«* -%-lf e 2a '{s + nh)™- 1 ds. 

2a 2,a A Jo 

From (71), h ^ 0, so that (93) yields 

./ <#(^ + />) e 2< "(* + n/O n/2 . (94) 
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Substituting (94) into (90), 



I 



•>, 



[«■ (± + *)(.+ a)] 



i/2 



ffi(-+i)(«) I < ^ Isf- + h) c- U S£2 -£- =L-. (95) 



© 



Recalling that » is some even integer ^ 2 in (95), (95) agrees with (75). 
Noting (79) and (82), the results of (75) and (76) hold for all n by 
induction. 

We now use the results of (75) together with (74) to obtain bounds 
on stability for those cases where the reflection coefficient r(z) is re- 
stricted as in (70). This analysis is almost identical to that of Section IV, 
(52)-(57), for the general case, modified by replacing the relation of 
(46) by that of (75). Thus, making the substitution 



I \r(s)\dsJ 



ni 






3- 1|; 



throughout (54)-(57), we obtain, corresponding to (57), the following 
sufficient condition for stability in the present case, after a minor modi- 
fication of the summation index : 



(£+*)£ 



#i± 



/,' (^ + /0 '-L, ' ml 



<--**#. (97) 



L N is the total length of the device. The summation of (97) is found in 
closed form by the analysis given in the Appendix. Using the final result 
oi the Appendix (137), the final results of this section may be sum- 
marized as follows: 

If the reflection coefficient r(z) (continuous, discrete, or a combination 
of both) satisfies the condition 

R-(z -/) ^f o \Hs)\d8gR-(z+g)i R>Q,f£0,g>0 (gga) 

h m f + g; h ^ 0. 
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then a sufficient condition for stability of the active line (with reflection) 



is 

where 



( 1+ 2^) ^ r 9r / u Sn\] 08b) 



(98c) 



\2a / c 

and n is given by 

fi = e iT \ n < e. (98d) 

The results of (98) are illustrated in Fig. 5, which shows the maximum 
value of R for which stability is guaranteed by (98) versus the nominal 
total gain 20 logio e aLrf , with 20 log 10 e ah as a parameter. 

A greatly simplified but slightly poorer version of the stability condi- 
tion of (98) may be obtained in the high-gain case. As one example, 
suppose the one-way gain of the active line exceeds 10 db, 

e *aL N ^ 1Qj g686 uLn £ 10db> aLs ^ 1151 . (99) 

If 8 satisfies the sufficient stability condition of (98b), it must also 
satisfy the weaker inequality 

S < *■*-«, (100) 

1 + 2atl 

Substituting (99) into (100), 

5 < 0.1. (101) 

From (98d), n is a monotonic increasing function of 8. Therefore 

n < 1.118. (102) 

Further, since from (98d) 

8n = ]nrt, (103) 

8ri is a monotonic increasing function of r x , so that 

8n < 0.1118. (104) 

Now writing out the right-hand side of (98b), 

exp \-2oLm (l + H)l - exp {-2aL N ) exp ("^ ^i) . ( 105 ) 
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0.10 



NOMINAL TOTAL GAIN, 20 LOG 10 e aL N IN DECIBELS 

Fig. 5 — Exact and approximate bounds on R for which stability is guaranteed. 



we investigate the exponent of the second factor on the right-hand side 
of (105). From (100), 



T 8ri< TT2ah 



N —2al.\ , ,-, i —2a Lid 

"•ri < 2aL N e "-r,. 



(100) 



The right-hand side of (100) is a monotonic decreasing function of 2aL N 
for 2aL N > 1. Therefore, substituting from (99) and (102), (100) yields 



^ 5r, < 0.2574. 



exp [-T 5ri ] 



> 0.7731. 



(107) 

(108) 
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Finally, using (104) and (108) in (98b), we obtain the following suf- 
ficient condition for stability, subject to (98a) ; 

R < 0.8287 - , 2 " t e~ aLN ; 8.686aL* ^ 10 db. (109) 
1 + 2ah 

The stability condition of (109) is slightly poorer than the stability con- 
dition of (98b), (98c), and (98d), from which it was derived. As the lower 
bound on the gain of the active line increases beyond 10 db and ap- 
proaches oo, the numerical factor 0.8287 in (109) increases and ap- 
proaches 1. Equation (109) or a similar result is useful in illustrating the 
general behavior; however calculations using the basic result of (98) are 
straightforward. Theresultof (109), with the numerical factor 0.8287 — > 1, 
is also shown as the dashed curves of Fig. 5, illustrating the way in 
which this approximate stability condition approaches the exact result 
of (98) in the high-gain case. 

VI. EXAMPLES AND DISCUSSION 

Consider first an active line with two discrete reflectors of equal mag- 
nitude c at the ends of the line, z = and z = L 2 . c is of course real; 
for convenience we assume c > 0. In this simple case the exact stability 
condition is readily found, and may be compared with the two bounds 
derived above. From (8) of Section I, the transmission gain of this de- 
vice in the stable region is 

G T = - , (110) 

x u 

where from (l)-(7) 

Xll _ e rL >(l + cV* rL '). (Ill) 

The condition for stability is readily found as described following (8) 
[this procedure is similar to that used in Section IV, (52)-(57), and 
Section V, (9G)-(97), in obtaining bounds on stability]. Replacing c 
by ec, where e is a numerical parameter greater than 0, and using (2), 

xu - e TL *[l + («e)V i-l V* tti ]. (112) 

For small enough t the minimum value of x u , and hence the maximum 
value of gain G T of (110), occurs at 

2/3L 2 = ±tt, ±3tt, ••• . (113) 
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Hence 

|. ril | min = e** L *[\ - ( e c) 2 e +2a/ ' 2 ]. (114) 

As e increases from zero, instability will take place at a value of e for 
which 

| X\\ |min = 0, 

, ,, a ,, 1 (115) 

(ec) c 2 = 1. 

Hence the original device (with e = 1 ) will be stable if 

c < e~ aLi . (116) 

Equation (116) is an exact condition for stability for the active line 
described above, with two equal reflectors at the ends. We now compare 
this exact result with the bounds described above. 

Consider first the bound of (9) or (62). This result is a sufficient con- 
dition for stability for any arbitrary distribution of discrete reflectors, 
and so must apply to the special case above. Setting N = 2, c\ = c% = c, 
this general bound guarantees stability if 

—aL 2 

tanrT 1 c < sinh -1 ^— . (117) 

V2 

Equation (117) yields 



I e ■ 

c < 



V2Vl + \c~' laLl 



(us; 



as a sufficient condition for stability for an active device with two equal 
reflectors of magnitude c at the ends. Comparing the bound of (118) 
with the exact stability condition of (116), we see that the general 
bound of (9) or (62) is conservative in the present special case; i.e., 
the device with two equal reflectors at the ends remains stable for the 
reflector magnitude c larger than that guaranteed by the general bound 
ol (9) or (62) by a numerical factor that varies from -y/s to y/2 as the 
gain aL 2 varies from to 00 . Therefore the general bound on stability 
given in (9) or (62) cannot be improved by a factor greater than \/2 
[i.e., this factor to multiply the right-hand side of (9) or (62)]; of course 
it may be that no improvement at all is possible, and that some distribu- 
tion of reflectors can be found for which (9) is satisfied as an equality at 
the boundary of instability. 

Next, consider the bound of Section V, (98), applied to the above 
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special case, i.e., two discrete reflectors of identical magnitude c at the 
ends of the active line. In (98) we set h = L 2 , R = (tanh -1 c)/L 2 , 
to yield the following (precise) bound on stability: 

8n < 2aU fl - tet , (n9a) 



where 



and r\ is given by 



1 - Sri 1 + 2aL 2 



5 = (tanh- 1 c) 2 - 1 +^ L2 (H9b) 



n = e 



or, 



n < c. (H9c) 



The bound on c for stability is readily determined numerically from 
( 1 19) as a function of aL 2 . However, when the one-way gain of the active 
line is large, aL 2 » 1, the bound of (98) takes on the form of (109), 
with the numerical factor 0.8287 — ► 1 since ah 2 » 1 (i.e., the gain is 
taken to be very large, not simply greater than 10 db). Thus the ap- 
proximate bound on stability in the present case becomes 

tanh" 1 c < : 2 °b T e- aL >; aL, » 1. (120) 

1 + l<xL 2 

The symbol < indicates that the relation of ( 120) is not a precise bound, 
but merely gives a good numerical approximation to the precise bound 
if aL 2 is large enough. Comparison of the (imprecise) bound of (120) 
with the exact stability condition of (116) shows that in the high-gain 
case, ah 2 » 1, where c « 1, the specialized bound of Section V, (98), 
yields bounds on the magnitude of the reflection c in the present special 
case (two equal reflectors at the ends of the active line) that approach 
those of the exact condition for stability. Consequently the bounds of 
(98) cannot be further improved (in their present form). 

The case of iV identical, equally spaced reflectors was studied in Sec- 
II of Ref . 1 , where simple expressions for stability were found in the high- 
gain case. If the total gain is large and the gain per section small, com- 
parison of (109) (with the factor 0.8287 -> 1) and (98a) with (43) of 
Ref. 1 shows again that the bound on stability of (98) cannot be 
further improved. It is of interest to see how close the bounds of (98) 
come to the exact value corresponding to instability in a few cases of 
interest. For this purpose we consider examples (i), (ii), and (in) of 



TRANSMISSION LINE STABILITY 



323 



Section II, Ref. 1. In (98) we set 

h = 1, R = 



tanh c 
l~ 



(121) 



and compute upper hounds on | c | that guarantee stability. It is also of 
interest to compare the general bound of (9) or (62) for this case. 
Table I summarizes these results. The bounds of (98) are quite good 
when the total gain is high, aL N ~S> 1, and when the gain corresponding 
to the distance I is small, ad « 1 ; for these conditions the stability condi- 
tion of (98) gives much better results than the more general stability 
condition of (9), because in the former we have made use of additional 
information regarding the distribution of reflectors. 

Table I — Identical, Equally Spaced Reflectors 

N = number of reflectors 

Gain (db) = 20 logio e Nal = 20 logio c l n = one-way gain of active line 
in db 

| c | mnx = maximum value of | c | for stability, as determined in Section 
II, Ref. 1 

Bound on | c j — (98) = maximum value of | c | for which stability is 
guaranteed by (98) 

Bound on | c | — (9) or (G2) = maximum value of | c | for which sta- 
bility is guaranteed by (9) or (62). 



Case 
(Sec. II, Ref. 1) 


N 


Gain, db 


C Imax 
(Sec. II, Ref. 1) 


Bound on | c | 
(98) 


Bound on | c \ 
(9) or (62) 


(«) 

(m) 


30 

300 

50 


30 

30 

5 


0.00860 

0.000860 
0.06501 


0.00590 

0.000710 

0.01105 


0.00149 

0.000149 
0.0130 



t Note that for this case in Ref. 1 the high-gain approximation given there was 
inappropriate, so that this result was obtained by use of a computer. 

Filially, we consider the application of the above stability conditions 
to some of the problems involving random reflectors studied in Ref. 1. 
The stability of the various deterministic cases discussed above in the 
present section has been treated exactly here or in Ref. 1 without using 
the new results of the present paper; these cases have been discussed in 
the present section both to show that any possible improvement in these 
general stability conditions must be quite small, and to provide partial 
confirmation of these results. However, the application of (9) and (98) 
to cases involving random reflectors provides the principal motivation 
for the present analysis, since no other information whatever is available 
regarding stability in these cases. 

Let us consider the example of the first part of Section IV, Ref. 1, in 
which the average normalized loss and the rms loss fluctuation were 
determined for an amplifier with reflections having identical magnitude 
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but random spacing. The following parameters were chosen for this 
illustration : 

l k = spacing between (k — l)th and fcth reflectors [(3) and Fig. 1] 

/ = (l k ) t average value of h , independent of k 

Ck = magnitude of A*th reflection coefficient [(3) and Fig. 1] 

c k = Co ; all reflectors identical, c > (122) 

N = 30, number of sections 

20 logio e""' = 30 db, nominal total gain 

20 logio e al ° = 1 db, nominal gain per section. 

The following assumptions were made in these calculations of Ref. 1: 

(a) U is always a large number of wavelengths; 

/3/a»2tt, /3/o»2tt. (123) 

(b) The distribution of the h about their mean l is very narrow with 
respect to the mean, but wide compared to 27r//3; further, this 
distribution is smooth and symmetrical about / . 

The probability density for h did not have to be further specified for 
the calculation of average loss and rms loss fluctuation in Ref. 1. (Note 
however that in the calculations of Ref. 1 for the covariance of the loss, 
the specific form of the probability density for h must be known, and 
was assumed to be Gaussian in Ref. 1.) The average loss and the rms 
loss fluctuation for the amplifier of (122) were given in Fig. 9 of Ref. 1 
versus c , the magnitude of the reflections. These curves were shown 
dotted for c > 0.00860, because it was known that instability is possible 
in this range, in particular for h = h , i.e., equally spaced reflectors [see 
Section II, Ref. 1 and case (»), Table I]. However it was noted that this 
was only a symbolic reminder of the unsolved question of stability; these 
results are valid for small enough c , but how small was not known from 
the results of Ref. 1 . 

We illustrate the utility of the results of the present paper by applying 
them to this problem; these results provide useful information concern- 
ing stability in this case, and of course in many similar problems. For 
convenience we make one further assumption in addition to those 
mentioned following (122): 

(c) The distribution of h about its mean h is strictly bounded; in 
particular 

\h - /o| ^ vk\ (124) 

further, we assume for convenience that 

v < 1. (125) 
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v is in (124) the upper bound on the fractional deviation in spacing from 
its average value; the restriction of (125) requires that h ^ 0, and so 
prevents the order of the reflectors from being altered. In practical cases 
we will be interested in small values of v, 

u«\. (126) 

We determine upper bounds on the reflector magnitude Co that guaran- 
tee stability, as a function of v, the maximum fractional deviation in 
spacing between reflectors. For v = the reflectors are equally spaced; 
Ref . 1 or Table I shows that stability is guaranteed if 

co < 0.00860, v = 0. (127) 

Next, the bound of (9) guarantees stability independently of the par- 
ticular distribution of reflectors. Since however the total length may 
vary somewhat, we must in (9) set 

L N m /, 3 o = 30fo(l + v), (128) 

yielding 

Co < 0.00149(0.03162) r (129) 

as a sufficient stability condition. 

Finally, we apply the bound of (98) to this example. We set 

B = tonp !) (130) 

«0 

h - (1 4- 60v)/ (131) 

and make use of (128) in (98) to obtain a sufficient stability condition. 
The sufficient stability conditions of (127), (129), and (98) are 
plotted in Fig. 6; the result of ( 129) is identified as originating from (9), 
and that of (127) from Section II of Ref. 1. The curves of Fig. 6 have 
been plotted out to fractional spacing variations v of 10 per cent; over 
this region the stability condition of (98) is superior to that of (9). 
However the bound of (9) [i.e., (129)] will be superior to that of (98) 
for large enough v. Note that the factor (0.03162)' in (129) arises from 
the fact that the total length and hence the total gain is subject to 
statistical fluctuation [a similar factor occurs in using (98) for the 
problem]; in the range of probable interest, i.e., for very small fractional 
spacing fluctuations v, this numerical factor will be close to 1. The fact 
that the limit of the bound of (98) as v — > is substantially below the 
maximum value of c given by (127) is due to the fact that the nominal 
gain per section in the example of ( 122) is 1 db, which is not too small; 
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as the gain per section decreases these two quantities will approach each 
other, as indicated above. 

These results, plotted on Fig. 6, show that the range of c over which 
the calculations of Section IV of Ref . 1 are guaranteed to be valid. If the 
maximum fractional variation in the spacing between reflectors is very 
small, then the results plotted on Fig. 9 of Ref. 1 are valid for c up to 
approximately 0.00590. 

The stability conditions of (9) and (98) may be applied to a variety 
of similar problems. In the above example we have found the maximum 
value of Co for which stability is guaranteed, i.e., for which the probability 
of oscillation is zero, as a function of the maximum departure of the 
spacing between reflectors from its average value. The results of (9) 
and (98) may also be used to determine an upper bound on the proba- 
bility of oscillation in similar problems where no absolute guarantee of 
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Fig. 6 — Bounds on magnitude of coupling coefficient to guarantee stability 
for amplifier of (122), with reflectors of identical magnitude and nominally equal 
spacing. 
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stability can be given, e.g., perhaps in cases where the probability dis- 
tribution for the spacing deviations is not strictly bounded. 

The main emphasis of the present paper has been on the discrete 
case; the continuous case was introduced only as an intermediate step 
leading to the desired results. However, it is clear that related problems 
with continuous reflection may be studied for stability using the general 
results derived above. 

Finally, the present calculations have assumed for definiteness a rather 
special model; i.e., the forward and backward gains have been assumed 
equal and a particular form has been taken for the matrix of the dis- 
crete reflectors. These assumptions are not essential to the analysis; 
similar results can be derived for many related cases of interest, such as 
systems using isolators to partially attenuate the backward waves, etc. 
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APPENDIX 

Summation of the Series S = zl ; 

«=o n! 

The summation of (97 ) was initially performed by a method suggested 
by S. 0. Rice, employing contour integration; this method is straight- 
forward but lengthy. A much shorter analysis presented by the unknown 
referee is given here. It has been shown that 6 

e" = 1 + t a(a - f r " y" (132) 

n=l Ti\ 

where 

y = xe bx and \yb\ < (1/e). (133) 

Differentiate (132) with respect to y and then set y = 1 to obtain 

(a— b)x oo r/_ 1 \ jin 

. _ [(a - 6) - m 

1 + bx ,,=o n\ 

x = e~ bz and |6| < (1/e). (135) 



where 
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Finally, set 



to obtain 



a = z - 5, b = -8, x = r, (136) 



[z + 8n) n e ri * Q 8 K 1 



fr'o n! 1 - an ' e 

where ?'i is given by (137) 

?'i = e iri , n < e. 
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